The rare decays of D mesons 

S. Fajfer'^''*, A. Prapotnik*^, S. Prelovsek'"'' , P. Singer'^, and J. Zupan^ 
'^Physics Department, University of Ljubljana, SI-1001 Ljubljana, Slovenia 
'^J. Stefan Institute, SI-1000 Ljubljana, Slovenia 

■^Department of Physics, Technion - Israel Institute of Technology, Haifa IL-32000, Israel 

The flavor changing transitions in the c — » wy, c — > ^77 and c — > ul'^l" ofi^er the possibility to search for 
new physics in the charm sector. We investigate dominant decay mechanisms in the radiative decays D V7, 
D — > P{V)l'^l~ , D — + 77 and we discuss chances to see physics beyond the standard model in these decays. In 
addition, we analyze Cabibbo allowed D Kii'y decays with nonresonant Ktx, and we probe the role of light 
vector mesons in these decays. 



The search for physics beyond the standard 
model has been focused on the down-like quark 
sector, while the up-like quark sector has been 
less researched. The rare D decays offer an op- 
portunity to investigate the FCNC effects in the 
charm sector. To date, no radiative or dilepton 
weak decay of D has been detected. Only upper 
bounds have been established so far for a sizable 
number of the radiative or dilepton weak decay of 
D mesons. The radiative decays p'^, uj'~', (p, 

-I-7 were recently bounded to branching 
ratios in the 10"'' range, which is approaching the 
standard model expectations (see, e.g. j|] where 
additional previous works are listed). The dilep- 
ton decays D Pl^l^, D are the sub- 
ject of intensive searches at CLEO and Fermilab 
m - Here again, with upper bounds of < 10~^ 
for branching ratios of the various modes one ap- 
proaches the expectations of the standard model 
§-§. The situation should improve in the future, 
due to new possibilities for observation of charm 
meson decays at BELLE, BABAR and Tevatron. 
Recently, upper limits in the 10~^ — 10~^ range 
were established also for dilepton decays 
with two nonresonant pseudoscalar mesons in the 
final state £1° (tt+tt- , K-Tr+ , K+K') n+ p' , 
though no comparable results are available yet for 
similar photonic decays. 

The inclusive c ^ wy transition is strongly 



*Talk given by S. Fajfer 



GIM suppressed at one loop electroweak order 
giving the branching ratio of the order 10~^^. In 
our approach we include the c — > U7 short dis- 
tance contribution by using the Lagrangian 

C^-AC^^^F^,[ua^^^\{l + ^,)c], (1) 

where rric is a charm quark mass and A — 
(Gf/V^) VusV;^. We follow |0| and we take 
= (-0.7-1- 2i) X 10-2. The branching ratio in- 
duced by this QCD corrected effective Lagrangian 
is BR(c — > uj) ~ 3 X 10^*. A variety of mod- 
els beyond the standard model were investigated 
and it was found that the gluino exchange dia- 
grams within general MSSM give the largest 
enhancement 

BR(c ^ m7)mssm _ (2) 
BR(c ^ uj)sM ~ ■ ^ ' 

Unfortunately, the long distance physics screens 
such effects since it usually dominates the decay 
amplitude The long distance contribution is 
induced by the effective nonleptonic |Ac| — 1 
weak Lagrangian 

£ = - Ay [ai?27^(l - 75)gigj-7M(l " 75)c 
+ 029^7^(1 -75)gj?i7^(l-75)c], (3) 

with Aij — {Gf/^/2) V*q.Vuqij accompanied by 
the emission of the virtual photon. Here qij de- 
note the d or s quark fields. The effective Wilson 



coefficients are oi = 1.2 and a2 = —0.5 1^. In 
our calculations of the long distance effects we use 
H,^,^ the theoretical framework of heavy meson 
chiral Lagrangian ||l^ . In the treatment of the D 
mesons rare decays we use the factorization ap- 
proximation for the calculation of weak transition 
elements. We consider the use of this approach 
to be justified by the "near" success of the ap- 
proach for the nonleptonic amplitudes. In Table 
1 we present the branching ratios of D ^ 
decays [||. The uncertainty is due to relative un- 
known phases of various contributions. Although 



enhanced up to 0(1) in the MSSM. In addition to 
the c — ^ U7 searches in the charm meson decays, 
we have suggested to search for this transition in 
i?*7 decay ||lj], where the long distance 
contribution is much smaller. 

The c ul~^l~ amplitude is given by the 7 
and Z penguin diagrams and W box diagram at 
one-loop electroweak order in the standard model, 
and is dominated by the light quark contributions 
in the loop. In Table 2 we give the c — > ul^l~ 
branching ratios calculated in the standard model 
and MSSM j^]. The amplitudes for exclusive de- 
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the branching ratios are dominated by the long 
distance contributions, the size of the short dis- 
tance contribution can be obtained from the dif- 
ference of the decay widths r{D'-* ^"7) and 
r(D° Lo^) jl^. Namely, the long distance 
mechanism cu dd-y overshadows the cu — > uu'y 
transition in D'^ p^^ and D'^ — s- ^7, the p'^ and 
uj mesons being mixture of uu and dd. However, 
the LD contributions are mostly canceled in the 
ratio 



R 



BR(i:>" pS) - BR(D0 W7) 



cx Re 



BR{D° W7) 
A{D° uu'y) 
A{D" dd-f) ' 



(4) 



which is proportional to 
A{D'^ uu'y) driven by c 
is 6± 15% in the standard model 



the SD amplitude 
-> U7. The ratio R 
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branching ratios. 
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cays D — + Vl'^l~ and D Pl'^l~ are domi- 
nated by the long distance contributions. The 
branching ratios for these decays, as obtained in 
1^,^, are given in Tables 3 and 4. In these Ta- 
bles the first column represents the SD contribu- 
tions, while in the second column the rates com- 
ing from the LD contributions are given. The 
rates for D Ve^ e~ are comparable to those 
in Table 3 and can be found in ||]. There is a 
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significant improvement of the experimental up- 
per bounds for D Pl^l^ decay rates, as re- 
cently obtained by FOCUS jj], with new upper 
bounds of 10"^ or less, close to the theoretical 
predictions [||J|]. The allowed kinematic region 
for the dilepton mass ma in the D Pl^l~ 
decay is mu = [2mi,mu — mp]. The long dis- 
tance contribution has resonant shape with poles 
at mil — nT'pO , ttIuj , m^f, [^,^ . There is no pole at 
mil = since the decay D P7 is forbidden. 
The short distance contribution is rather flat. 
The spectra of I? — + Pe+e~ and D Pfx^ ix~ 
decays in terms of mu are practically identical. 
The difference in their rates due to the kinematic 
region mu = [2me-, 2m^] is small and we do not 
consider them separately. 

The differential distribution for _D+'° 
7r+°Z+Z~ 1^,^ indicates that the high mass dilep- 
ton region might give an opportunity for de- 
tecting c — > ul^l^ . Before making a definite 
statement on such possibility, we should examine 
this kinematical region of high dilepton mass in 
D ttI^I^ decays more closely. For instance, in 
this region the excited states of the vector mesons 
p, u> and 4> may become important. Even with the 
presence of these states we found that the only 
viable channel for investigating the c — s- ul^l^ 
transition is D — s- nl'^l", in the kinematic re- 
gion of the dilepton invariant mass mu above the 
resonance (f>^ where the long distance contribu- 
tion is reduced. The kinematics of the processes 
D — + Vl^l^ would be more favorable to probe the 
possible supersymmetric enhancement at small 



mil, but the long distance contributions in these 
channels are even more disturbing |6|. 

Motivated by the experimental efforts to ob- 
serve rare D meson decays [p^ , and noticing that 
Bs 77 offers possibility to observe physics 
beyond the SM, we undertook an investigation 
of the Z?" 77 decay ||l^. The short dis- 
tance contribution is expected to be rather small, 
as already encountered in the one photon de- 
cays [0 , hence the main contribution would come 
from long distance interactions. The total ampli- 
tude is dominated by terms proportional to ai 
that contribute only through loops with Gold- 
stone bosons. Loop contributions proportional 
to 02 vanish at this order. We remark that the 
contribution of the order 0{p) does not exist in 
the 77 decay, and the amplitude starts 

with contribution of the order 0{p^). The chi- 
ral loops of order 0{p^) are finite, as they are 
in the similar case of K 77 decays. At this 
order the amplitude receives also an annihilation 
type contribution proportional to the 02 Wilson 
coefficient, given by the Wess-Zumino anomalous 
term coupling light pseudoscalars to two photons. 
Terms which contain the anomalous electromag- 
netic coupling of the heavy quark Lagrangian are 
suppressed compared to the leading loop effects 
fl^ ] . The invariant amplitude for — s- 77 decay 
can be written using gauge and Lorentz invari- 
ance in the following form, 



M = 



^ ki-k2' 
+M(+)e^''"^fci„A:2/5lei^e2., 



(5) 



where M^~'> is a parity violating and M^^^ a par- 
ity conserving part of the amplitude, while fci(2), 
^1(2) a-i'e respectively the four momenta and the 
polarization vectors of the outgoing photons. We 
give in Table 5 the numerical results for the am- 
plitudes originating from the different contribu- 
tions. Within this framework, the leading con- 
tributions are found to arise from the charged tt 
and K mesons running in the chiral loops and 
our calculation predicts that the D — > 27 decay is 
mostly a parity violating transition |l^ . We esti- 
mate that the total uncertainty is not larger than 
50% 1 16 , including possible effects from light res- 
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onances like p, K* , ao(980), /o(975). Accord- 
ingly, we conclude that the predicted branching 
ratio is 

BR(L>° 77) = (1.0±0.5) X 10^^ (6) 

We look forward to experimental attempts of 
detecting this decay. Our result suggests that the 
observation of D ^ 27 at a rate which is an order 
of magnitude larger than ( |q), c ould be a signal for 
the type of "new physics" [^,0 . 

In addition to the D meson rare decays in which 
the FCNC transitions might occur, we undertook 
a study of the Cabibbo allowed radiative decays 
A'°7r+7 and Z?" K^tt^'j with nonreso- 
nant K tt, which we consider to be the most likely 
candidates for early detection. Here we used the 
heavy quark chiral Lagrangian supplemented by 
light vector mesons, as the theoretical framework. 

These decays are the charm sector counterpart 
of the K — !■ 7r7r7 decays, which have provided a 
wealth of information on meson dynamics. Using 
the factorization approximation for the calcula- 
tion of weak transition elements we use the infor- 
mation obtained from semileptonic decays 0,|8| . 
The nonleptonic D —^ Kit amplitude cannot be 
calculated accurately in the factorization approxi- 
mation from the diagrams provided by our model. 



Such a calculation gives a rather good result for 
the Z?+ K'^tt^ channel but is less successful for 
the _D" K^TT^ decay. In order to overcome this 
deficiency and to be able to present accurately the 
bremsstrahlung component of the radiative tran- 
sition, we shall use an alternative approach for 
its derivation. It means, we take the experimen- 
tal values for the D Ktt amplitudes, in the 
calculation of the bremsstrahlung component. In 
order to accommodate this, we write the decay 
amplitude as 

+ F,[{q ■ e){p ■ k) - [p ■ e){q ■ k)] 

+ F2e^''^''e^v^k0q^}, (7) 

where Fq is the experimentally determined D 
Ktt amplitude and Fi, F2 are the form factors 
of the electric and magnetic direct transitions 
which we calculate with our model. When in- 
termediate states appear to be on the mass shell, 
we use Breit Wigner formula. Thus, we get for 
the branching ratios of the electric transitions in- 
cluding bremsstrahlung, with |i^o| determined ex- 
perimentally and taking the photon energy cut 
> 100 MeV 

BRiD+ ^ K\+j)§^^^^ 

= (2.3-2.5) X lO"''. (8) 

For the D'^ radiative decay we get 

BR{D" ^ /^"^+7)pv,ex 

= (4.3-6.0) X lO"''. (9) 

The uncertainty in the Fq/Fi phase is less of a 
problem in £)+ K'^tt^j than in D'^ — K~Tr~^'y. 
If we take the bremsstrahlung amplitude alone 
as determined from the knowledge of |Fo|, dis- 
regarding the direct electric Fi term, the above 
numbers are replaced by 2.3 x 10""* for _D+ decay 
and 5.5 x 10^"* for the Z?" decay. The contribu- 
tion of the direct parity violating term (putting 
Fo = 0), is BR(i:i+ X°7r+7)dir,pv = 1.0 x 10"^ 
and BR{D'^ /^"7r+7)dir,pv = 1.64 x lO^^. For 
the parity conserving direct magnetic transition 



we get 
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BR{D° 
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-5 



(10) 

(11) 



Hence, the two direct transitions are predicted to 
be of about the same strength. 

If we disregard the contribution of vector 
mesons to the direct part of the radiative de- 
cays, the parity-conserving part of the amplitude 
is considerably decreased, by one order of magni- 
tude in the rate in Z3+ K'^tt'^'j decay and by 
two orders of magnitude in —f K^n^j. On 
the other hand, their contribution is not felt in 
a significant way in the parity - violating part of 
the amplitudes. In any case, the detection of the 
direct part of these decays at the predicted rates, 
will constitute a proof of the important role of the 
light vector mesons. 

We conclude by expressing the hope that the 
interesting features which these decays provide 
and were analyzed in JlOt , will bring to an ex- 
perimental search in the near future. 
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